Hubble diagram of gamma-ray bursts: Robust evidence for a Chaplygin gas 
expansion-driven universe with phase transition at z ~ 3 
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The Hubble diagram (HD) of Gamma-Ray Bursts (GRBs) having properly estimated redshifts is 
compared with the predicted one for the Chaplygin gas (CG), a dark energy candidate. The CG 
cosmology and that of Friedmann and A-CDM models are studied and confronted to the GRBs 
observations. The model-to-sample x 2 statistical analysis indicates the CG model as the best fit. 
The present GRBs HD plot exhibits a marked trend: as one goes back in time, it gets much closer 
to the predict HD for a Friedmann universe. This clear trend conclusively demonstrates that a 
transition from decelerate to accelerate expansion did take place. However, contrarily to claims 
based on supernovae type la, the transition redshift lies somewhere between ~ 2.5 < z ~ 3.5 rather 
than at z ~ 0.5 — 1. All of these striking features of the GRBs HD constitutes the most robust 
demonstration that the Chaplygin gas can in fact be the universe's driving dark energy field. 
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Introduction. — Observations of supernovae type la 
(SNIa) have led to the current view that our universe un- 
derwent a late-time transition to accelerate expansion at 
a redshift z ~ 1. The driver of such unexpected dynamics 
is an exotic component of the universe's content dubbed 
dark energy (DE), a smoothly sparsed energy field with 
no familiar counterpart among the currently known forms 
of matter-energy. Lots of theories have been conceived to 
explain this striking phenomenon in contemporary cos- 
mology. One of these candidates for DE is the so-called 
Chaplygin gas (CG), a strange fluid described by a quite 
unsual equation of state (EoS) 

P=-~, (1) 
P 

where p represents the pressure, p the fluid density and 
A is a constant. It came after the russian aerodynamicist 
Chaplygin who in 1904 brought it in to explain the lifting 
force on a plane wing in some aerodymamic phenomena 
0. By 2001, Kamenschik, Moschella and Pasquier Q 
recognised its relevance to cosmological studies, in par- 
ticular, with respect to the claimed cosmic acceleration. 
They showed that the CG model exhibits excelent agree- 
ment with observations. Besides, the model predicts a 
larger value for the effective cosmological constant. The 
same researchers noticed then that the model can be gen- 
eralized in the form: p = — -4r, and consider the case with 
the density power-law exponent a = 1/3 |2j . It was then 
realized that the CG EoS has a clearly stated connection 
with string and brane theories[3j,[7|,[8]. 30] 

Here we build the Hubble diagram (HD) of gamma-ray 
bursts (GRBs)j3l|, which already reach redshift z ~ 7, 
and SNIa, and confront it with predictions of the cosmo- 
logical model in which our universe is filled-in with this 
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FIG. 1: (color on-line) HD of GRBs and the CG model of 
DE. The plot demonstrates that a phase transition did take 
place at a redshift 2.5 < z ~ 3.5, or even earlier (right panel) 



sort of DE, the so-called Chaplygin gas. In comparing 
the CG predicted HD with both GRBs and SNIa data is 
verified that the best fit to the observational data clearly 
correspond to this exotic fluid, when compared to the HD 
for the Friedmann-Lemaitre-Robertson- Walker (FLRW, 
the universe contents only matter) and the Lambda Cold 
Dark Matter (A-CDM, Friedmann cosmology with A) 
scenarios. Three major achievements are attained in this 
investigation: a) for the first time is presented the HD 
of GRBs in confrontation to the CG theoretical predic- 
tion, and that of the FLRW and A-CDM models (see 
FIG. 1). b) it is shown that the best fit (x 2 -statistics) 
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to the GRBs data is provided by the Chaplygin gas sce- 
nario, c) the resulting HD clearly exhibits a transition 
from a Friedmann-dominated to a late-time accelerating 
universe, with the transition taking place at a redhsift 
around 2.5 < z ~ 3.5 (see FIG. 1), and driven by the 
CG. Besides, the similar analysis for SNIa allows one 
to verify that the HD for SNIa clearly violates the CG, 
FLRW and A-CDM expansion law. This disagreement 
with the Chaplygin gas HD, once again suggests that 
perhaps there is something wrong with the SNIa obser- 
vations and that its astrophysics deserves to be revisited. 

Chaplygin gas cosmology. — The enthusiasm about the 
Chaplygin gas dynamics stems from the possibility of uni- 
fying dark matter and dark energy in a cosmic fluid that 
could be represented by a scalar field. As one can notice 
from the evolution law for its energy density 0, 

where a is the scale factor of the universe and B is an 
integration constant, the CG corresponds to the A-CDM 
scenario for the parameter a = 0, and A s = A/ p)£ a = 
1. A direct analysis of Eq.(J2| shows that the present 
acceleration phase, which should have taken place at a 
redshift 



z t = [3u(zt) + 1] 
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n, 
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(3) 



where u>(z a ) 



Jo *\+J dz , leads to an asymptotic 



ln(l+z) JO l+z' 

(a — ► oo) stage where the EoS is dominated by a cosmo- 
logical constant (SnA 1 /^ ), whereas at earlier epochs 
the energy density is dominated by a non-relativistic mat- 
ter. Notice also that the value at of the scale factor that 
signalizes the start of the late-time acceleration phase is 
given by the roots of the equation [25| 



d = a(H + H 2 ) , 
and is related to the redshift value z t such that — 



(4) 



l 

a a l+z t ' 

The highlighted dual behavior is the basis of the unifi- 
cation scheme provided by the CG model. This unifica- 
tion becomes possible if one takes benefit of a complex 
scalar field $ of mass that admits an inhomogeneous 
generalization, as demanded by models of structure for- 
mation in the universe, and is described by a Lagrangian 
density 0, @ 

C = gva ,*»V-V(\*\f (5) 

where $ = (tj)/\/2 m) exp(— ira^ff). The cosmic dynamics 
of this field was developed in Ref.0,0]. 

As pointed out above, the GCG Model can be idealized 
by a perfect fluid with an EoS given by 

P=~— , (6) 

where A and a are constants. When a — 1 we re-obtain 
the EoS for the CG scenario (Eq^J. In principle, the 



parameter a is restricted in such a way that < a < 1. 
However, possible values for a ^ (0, 1] are considered in 
the accompanying paper 9|, where we also analyze the 
effects of imposing the energy conditions to the cosmic 
dynamics of the CG as to be compared with GRBs and 
SNIa observations. 

The universe content can be envisioned as having a 
pressureless matter, needed to account for the presence 
of baryons in it, and also dark matter (also pressureless) 
and dark energy making-up the CG. Hence, the dynamics 
of the Universe is worked out through the Friedmann's 
equation and the evolution equations for non-interacting 
baryonic matter and Chaplygin gas[25j 
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(8) 
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where p m and p c stand for the pressureless matter and 
Chaplygin gas component, respectively. As usual, k = 
0, 1, —1 indicates a flat, closed and open spatial section. 
The conservation law for each of these fluids H8I9|) 



reads: 



and p c = (A+ a3{1+a) ) 



l/(l+a) 



respec- 



tively. The value of the scale factor today is taken equal 
to unity, ao — 1. Hence, p m o and p C Q — (A+B) 1 ^ 1+a ^ are 
the pressureless matter and GCG densities today. Elimi- 
nating from the last relation the parameter B, the GCG 
density at any time can be re-expressed as 



Pc = PcO [A + 
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(10) 



where A — A/ p C Q. This parameter A is connected with 
the sound velocity for the Chaplygin gas today by the 
relation |^ = v 2 s = aA. 

Our main purpose here is the theoretical distance mod- 
ulus vs. redshift relation, i.e., the HD of the CG model, 
to compare it with GRBs observations. For this, we need 
the luminosity distance |2lL l2^| 



A °0 
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(11) 



with n the co-moving coordinate of the source. As light 
propagates on a null geodesic, i.e., 



ds 2 = c 2 dt 2 



2 dr 2 



l-kr 2 







(12) 



the Friedmann's equation Q allows to re-cast the lumi- 
nosity distance as 



d L = (l + z)S[f(z)} 



(13) 



where S(x) = x (k — 0), S(x) = sinx (k — 1), S(x) = 
srnhx (k = — 1) , and the function f(z) being given by 
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FIG. 2: (color on-line) HD of GRBs and the CG model of DE for the 52 (first pair of plots) and 24 GRBs (second pair) 
samples, including the case of a cosmological constant with Qa = 0.73. Also the residual HD of the CG vs. Friedmann models 
confronting the 52 GRBs sample, and the HD small variation for several values of the CG a-parameter are presented. 
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with the definitions 



£2 , ■ 



8ttG p m Q 



8ttG Pcq 



Hi 



n 



feO 



Hi 



(15) 



such that the condition Q m o + ^co + Qka = 1 holds. The 
final equations did use of the redshift vs. scale factor 
relation: 1 + z = ^ a °~- >1 . 

Data analysis and results. — For the present analysis 
we have taken benefit of two samples of GRBs as com- 
piled by Schaefer 0], 52 GRBs that were analyzed by 
five different methods, and by Bloom, Frail and Kulka- 
rni ^3]> 24 GRBs. Both having properly estimated dis- 
tance modulus, n(z), and redshifts. We also used the 
SNIa data as collected in the GOLD s amp le |24j. and 
the SuperNova Legacy Survey (SNLS) [23j. Our main 
results are collected in Figures 1, 2 and 3 (color on-line). 
They were obtained upon integration of the luminosity 
distance function in Eas. 1)13114(1 . for a flat universe with 
!l m o = 0.05 and £l c o = 0.95, and for a cosmological 



constant-like CG with fl c o = 0.73 and £! m o = 0.27. 

Conclusions. — We presented the Hubble diagram of 
a large sample of GRBs having properly estimated red- 
shifts together with the equivalent plot predicted by the 
cosmology of the Chaplygin gas, one of the various candi- 
dates to make-up the dark energy; the exotic fluid driving 
the universe late-time accelerated expansion, and FLRW 
and A-CDM models. After doing the cosmology of this 
model, we performed the statistical analysis and com- 
puted the model-to-sample x 2 . It was found that the 
CG definitely fits much better the HD of GRBs used in 
the present study, in comparison to the Friedmann and 
the pure A-CDM models. The present GRBs HD plot ex- 
hibits a marked trend indicating that as one goes back in 
time, i.e., to very high redshifts, it gets much closer to the 
predict HD for a Friedmann universe. This clear conver- 
gence demonstrates that a transition from decelerate to 
accelerate expansion did take place. However, and con- 
trarily to claims based on supernovae type la (SNIa), the 
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observations do not face. 



FIG. 3: (color on-line) Hubble diagram of SNIa: GOLD sam- 
ple (leftish pair) and LEGACY sample (rightish pair), and 
the prediction of the CG model of DE, including the case of a 
cosmological constant with Qa = 0.73. Similar plots are also 
given for a = 0.1, 0.3, 0.5, 0.7 . 



transition redshift lies somewhere between 2.5 < z ~ 3.5 
rather than at z ~ 0.5 — 1. Therefore, this GRBs HD con- 
stitutes the most robust demonstration that the Chap- 
lygin gas can indeed be the universe's driving dark en- 
ergy field. "En passin" , as still the GOLD and LEGACY 
SNIa HD locate far-above the CG predicted expansion 
law, this is further evidence of their intrinsic problems: 
or the Phillips relation should be revisited |2j], or per- 
haps there is some interrelation between the estimated 
luminosity distance residuals and internal extinction of 
the host galaxy, as some researchers have pointed out 
recently |28j . Finally, it is worth to quote that a sim- 
ilar conclusion regarding the transition redshift z ~ 3 
was achieved quite recently by Amendola, Gasperini and 
Piazza 26J upon the admission that dark matter and dark 
energy can strongly interact and evolve through a scaling 
regime pdm ~ Pde ~ a~ 3 ( 1+u,etf ' , with w e a a constant. 
However, although the present analysis here do confirm 
the transition around z ~ 3, as we pointed out here, 
and in Refe-HEl, the GOLD and LEGACY SNIa ob- 
servations seem to violate the general relativistic energy 
conditions (see Figure 3) for the (generalized) Chaply- 
gin gas and Friedmann models, a problem that GRBs 
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gas (GCG) model with the a-parameter since it will be 
the base of a detailed study that will be presented in 
the accompanying paper to this Letter|i|- Nonetheless, 
in the discussion ahead we will restrict ourselves to the 
case wherein u = 1. 

[31] Gamma-ray bursts (GRBs) are the biggest explosions in 
the universe. The 1997 february 28 Beppo-SAX discov- 
ery of the X-ray afterglow of a gamma-ray burst (GRB) 
allowed the first precise determination of the redshift of 
a GRB. This major breakthrough came to confirm the 
long-standing suspicion waving in the high energy as- 
trophysics community that GRBs arrive from, and their 
sources lie at, cosmological whereabouts. Since then, the 
possibility of using them as actual cosmological probes 
has estimulated the search for self-consistent methods of 
bringing GRBs into the realm of cosmology. A handful 
of attempts have been on trial. Since the introduction 
of the Amati relation ll], other largely promising tech- 
niques have appeared in an attempt to turn GRBs into 



reliable cosmological probes. These include the Ghirlanda 
relation the Liang-Zhang relation and the re- 
cently discovered Firmani et al. relation [14| . all of which 
taking into account the most relevant physical proper- 
ties of GRBs as the peak energy, jet openning angle, 
and both time lag and variability. Such discoveries hints 
at the long-sought Holy Grail of creating a cosmic ruler 
from GRBs observables to be achievable. Presently, after 
the first series of controversial statements on the viabil- 
ity of granting to GRBs the status of standard cosmic 
rulers |15l 1171 Hs|l . a definite consensus appears to be 
arising, and the hope to have a Holy Grail to do cosmol- 
ogy upon GRBs is renascenting [12, Ha, Ea 13 ■ Besides, 
we bring to the reader's attention a historical fact that 
has a clear correspondance with the present situation re- 
garding the GRBs cosmology. Despite having still large 
error bars, and in several cases a not so clear estimate of 
the redshift of some events, the present state of the cos- 
mology based on GRBs emulates the days during which 
Hubble discovered the expansion of the universe by using 
inhomogeneous and badly calibrated data from the neb- 
ulae he, Humason and others had observed IHJ. Nobody 
nowadays thinks of that his analysis as null and void, in 
spite of his methods having been not so standardized. 
Keeping these arguments in mind one can be confident 
of the worthiness for cosmological studies of the analysis 
being presented here. 



